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Abstract 

We have developed novel instrumentation for making Lagrangian measurements of temperature 
in diverse fluid flows. A small neutrally buoyant capsule is equipped with on-board electronics 
which measure temperature and transmit the data via a wireless radio frequency link to a desktop 
computer. The device has 80 dB dynamic range, resolving milli-Kelvin changes in temperature 
with up to 100 ms sampling time. The capabilities of these "smart particles" are demonstrated 
in turbulent thermal convection in water. We measure temperature variations as the particle is 
advected by the convective motion, and analyse its statistics. Additional use of cameras allow us 
to track the particle position and to report here the first direct measurement of Lagrangian heat 
flux transfer in Rayleigh-Benard convection. The device shows promise for opening new research 
in a broad variety of fluid systems. 

PACS numbers: 47.80.-v (Instrumentation for fluid flows); 44.27.+g (Convective heat transfer) 



1 



I. INTRODUCTION 



Scalar mixing in turbulent flows plays a crucial role in uncountable natural, medical, 
and industrial systems: spread of pollutants by wind and water; oceanic and atmospheric 
thermal convection; the life cycles of plankton; mixing in combustion engines and chemical 
reactors. A natural approach to understanding these examples is to measure the trajectories 
of the pollutant molecules, warm fluid elements, planktonic organisms, and chemical species, 
respectively, as well as the properties of the flow along these trajectories. This Lagrangian 
approach to turbulent mixing has been advanced signiflcantly with numerical and theoretical 
models l|, but experimental works are rare due to the difficulties in (1) distinguishing the 
identity of a fluid particle along its trajectory, and (2) making measurements along these 
paths . A few recent experimental studies have successfully performed resolved measure- 
ments of the trajectories of small solid tracer particles in turbulence 
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. Atmospheric 



research groups 



6| and oceanographers routinely employ meter-sized Lagrangian floats to 



study large scale flows Lagrangian probes have a distinct advantage over Eulerian probes 
which record at flxed positions in the flow: the measurements provide information about 
local physical processes as experienced by the fluid particle. 

With the aim of measuring Lagrangian scalar quantities in well-controlled laboratory 
flows, we have developed a miniature, wireless, neutrally-buoyant instrument, dubbed a 
smart particle. We present temperature measurements obtained with this device in turbu- 
lent thermal convection. Similar instruments have been used in medical applications Q], 
but to our knowledge we report the flrst such measurement device in the context of fluid 
dynamics. Simultaneous measurements of temperature and position have been made before 
in convection experiments, most recently by tracking optically the motion of microspheres 
containing thermochromic liquid crystal (TLC) j^. The shape evolution of thermal plumes 
very close to the thermal and hydrodynamic boundary layers was quantifled, but this method 
is not capable of tracking individual particles over large distances. Our smart particle mea- 
surements complement those of TLC particles by vastly improving temperature and time 
resolution as well as allowing observation of very long particle trajectories, though the smart 
particles are too large to investigate the flow in the boundary layers. 
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II. DEVICE DETAILS 



The details of the instrument are describe in the following five subsections. First we 
present a block diagram and general description of the whole system. Then, the radio 
frequency (RF) data transfer is described in detail. Next, we focus on the power management 
system. We then discuss the measurement and processing of temperature data and end with 
details about the position measurement technique. 

A. Overview 

The smart particle consists of a D = 21 mm diameter capsule containing temperature 
instrumentation, an RF emitter, a battery, and an on-off switch - illustrated in Figs. l(a,b). 
A resistance controlled oscillator is used to create a square wave whose frequency depends 
on the temperature of several thermistors. This square wave is used directly to modulate 
the amplitude of the radio wave generated by the RF emitter. The entire mobile circuit is 
powered with a coin cell battery and may be put in a low power standby mode using an 
externally applied magnetic field. 

The stationary parts of the system include an antenna, an RF receiver, 2 RF amplifiers, 
a high speed data acquisition system, and a desktop computer running Labview (Fig. la). 
The broad-band radio frequency amplifiers increase the voltage amplitude of the signal by 
a factor of 26 dB. The receiver is carefully tuned to demodulate the signal produced by the 
emitter. The receiver outputs a square wave identical to that generated by the resistance 
controlled oscillator. The frequency of this square wave, and hence temperature, is recovered 
on-the-fly using a Labview algorithm. In addition, the particle trajectory is recorded with a 
digital video camera, resulting in synchronous measurements of the position and temperature 
of the particle as it is carried about by the fluid. 

B. RF data link 

The radio frequency emitter is a MAX 7044 (from Maxim Integrated Products), which 
employs on-off keying (OOK) amplitude modulation. The RF carrier frequency is 315 MHz. 
The MAX7044 expects a digital modulating signal with a baud rate from to 100 kHz. We 
implement an 8 mm diameter, 8-turn coil emitter antenna with about 200 nH inductance 
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and 4 Q resistance. We use a split-capacitor impedance matching network to maximize 
the efficiency since the MAX7044 is optimized for a 125 Q antenna. The emitter draws an 
average of 5 mA and requires 2.1-3.6 V from the battery in normal operation. The emitter 
consumes 10 to 100 times more power than any of the other on-board components. 

On the receiving end we use the MAX1473 superheterodyne receiver, which demodulates 
OOK data and is designed to work with the MAX7044. The receiving antenna is a 315 MHz 
1/4 wave whip (Linx Technologies, Inc., ANT-315-CW-HD). Between the antenna and the 
receiver are two low noise amplifiers each with a gain of 13 dB at 315 MHz (MAX2640). 

C. Power management 

The size of the instrumentation capsule is minimized so that it capable of probing as 
small a spatial scale as possible. The component of the system which limits the size of 
the capsule most is the battery. It is a CR1616 lithium coin cell (Panasonic), which is 16 
mm in diameter and 1.6 mm thick. This battery can supply the necessary power for about 
3 hours. The circuit is also equipped with a magnetic field triggered switch so that the 
instrumentation may be turned off when we are not ready to acquire data. The switch is 
composed of a TLE4913 low power hall switch (Infineon Technologies) and a NL17SZ74 
single D flip-flop (ON Seminconductor). The hall switch outputs a high logic level when in 
the presence of a sufficiently large magnetic field. Feeding this output into the CP pin of the 
ffip-flop while D is connected to Q we have a "push-button" switch, which is operable with a 
hand-held permanent magnet from a distance of several centimeters. When the switch is in 
the "off' position, the circuit consumes less than 100 fiA. The hall switch and the fiip-fiop 
respectively require 2.4-5.5 V and 1.65-5.5 V from the battery. 

D. Temperature detection 

The system measures the spatially averaged temperature around the smart particle 
body. This is accomplished using four thermistors to set the frequency of a LMC555 
timer /oscillator (National Semiconductor) operating in astable mode. The thermistors pro- 
trude from the capsule wall by about 0.5 mm and are spaced evenly around its middle. 
The LMC555 outputs a square wave at logic levels with a period 1// = (-Ri + R2 + 2i?3 + 
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2i?4)C'/l-44, where Ri are the resistances of the four thermistors, and C is 47 pF. We employ 
0.8 mm, 230 kf2 thermistors (Epcos B57540G0234) with a response time of about 0.06 s in 
water, which is faster than the fastest temperature time scales of the flow we study. In the 
range of temperature 26 — 34 °C, the frequency of the square wave is in the range 22 — 26 
kHz; the sensitivity is 513 Hz/°C. The relationship between temperature and frequency is 
linear within 2 percent. The LMC555 requires about 100 fiA at a supply voltage between 
1.5 and 12 V. 

The RF receiver recovers the square wave signal, which is then recorded with a high-speed, 
14 bit, analog-to-digital converter (ADC) (National Instruments 5621 in a PXI system). 
Running on a desktop PC, Labview is used to control the data acquisition. The frequency 
of the demodulated square wave is computed on-the-fiy using a standard Labview library. 
Each measurement of frequency is computed from 100 ms of the square wave, i.e. about 
3000 periods, sampled by the ADC at 10 MHz. The resulting measurement resolution is 
about ±2 Hz, which is ±4 mK when converted to temperature. The sampling rate is about 
10 Hz. With maximum flow velocities in the range 1-2 cm/s and a particle size of 21 mm, 
we are oversampling the dynamics by a factor of order 10. 

E. Position measurement 

In order to measure the smart particle position, we use a standard webcam, interfaced 
with a desktop computer. A uniform well-illuminated image background is achieved with a 
back-lit sheet of frosted glass. We use Matlab scripts to control the camera and process the 
video data. First, each video frame is converted into a 2D binary array, using an adaptive 
threshold. Then, the particle position is extracted using an image recognition algorithm. 
The effective sample frequency is around 5 Hz. The resolution of the camera is 640 x 480 
pixels, so that the particle position is determined with 0.1 mm precision. 

III. MEASUREMENTS IN RAYLEIGH-BENARD CONVECTION 
A. Convection apparatus 

Our experimental setup, shown in figure Ic, is a traditional rectangular Rayleigh-Benard 
cell with height H = 40 cm and section 40 cm x 10 cm. The fluid is water and the walls of 
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the vessel are 25 mm thick PMMA. The upper plate of the cell is temperature regulated by 
a chilled water bath. The bottom plate is heated by 5 resistors, regularly spaced. Complete 
experimental details can be found in 101 ]. 

In the results reported here, the power input is P = 230 W, and the upper plate is 
regulated to Tup = 19°C, corresponding to a temperature difference between the top and 
bottom plates of AT = 20.3°C. As a result, the Rayleigh number is 



gpATH'^ 

UK 



Ra = — = 3.07 X 10^°, (1) 



where g is acceleration due to gravity, (3 = 2.95 x 10~^ is the thermal expansion 
coefficient of water and u = 8.17 x 10~^m^s~^, k = 1.48 x 10~''m^s~^ its viscosity and 
thermal diffusivity (values are given for a mean temperature equal to 29.1°C). The Nusselt 
number, measured as the total heat flux normalized by kATJL, is Nu = 167.9 ±0.2. Under 
these conditions, the convective regime is fully turbulent 10|, ll3| and the mean flow is a 
steady, system-sized, single convection roll with a rotation period of about 100 s. The 
characteristic thickness of the thermal boundary layer is ix ~ ^HNu^^ ~ 1.2 mm and that 
of the hydrodynamic boundary layer is £u ~ ixii^/ f^)^^'^ ~ 2 mm. Thus, the particle is too 
large to penetrate the boundary layers. 



B. Temperature measurement 

The particle and fluid density are carefully matched. This is achieved by initially adjusting 
the particle mass to within 1 percent of pfiiD'^/G, where pf is the density of water. Then, to 
obtain neutral buoyancy within 0.05 percent, p/ is finely adjusted by the addition of small 
amounts (about 1000 ppm) of pure glycerol, which has a density 20 percent greater than 
water. The particle is inevitably slightly lighter than the cool fluid near the upper plate 
and slightly heavier than the warm fluid near the bottom plate. Nonetheless, the particle 
explores all regions of the vessel without a clear bias as can be seen in figure 2b. 

We show in figure 3a, a time series of the temperature 6{t) recorded by the particle, 
during the trajectory shown in figures 2b, c. Apparent in the time series is a nearly periodic, 
large amplitude fluctuation, which reflects the particle's entrainment in the mean flow. This 
periodicity, as well as the smaller amplitude turbulent fluctuations are made more clear in 
the power spectrum shown in figure 3b. The spectrum also reveals a range of time scales, 
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corresponding to frequencies 0.016 < / < 0.230 Hz, which is consistent with a self-similar 
9"^ ~ /~" scaling, with an exponent a close to -2. This feature is in sharp contrast with the 
Eulerian temperature spectrum obtained when the particle is constrained to remain at fixed 
location in space. In that case, we have verified that the spectrum shows an /"^''^ scaling 
range. This value is in agreement with other experimental Eulerian studies of convection 
(e.g. provided the frequency spectrum be interpreted as a spectrum in space, which is 

plausible considering the persistent mean flow (i.e. Taylor's hypothesis). Bolgiano-Obhukov 
scaling also leads to a fc"''/^ wavenumber spectrum for temperature variations in real space. 
By contrast, the spectral behavior found here for the moving particle is reminiscent of 
velocity data for Lagrangian tracers in turbulent flows. Indeed when fluid turbulence is fully 
developed, a well-established feature of Lagrangian dynamics and Kolmogorov's theory is 
that the Lagrangian velocity spectrum should have an inertial range [s], [l^ - as we 
have verified also verified here from our independent measurement of the particle position. 
The upper cut-off frequency apparent in figure 3b around (~ 0.3 Hz) corresponds to the 
characteristic time of motion of the particle across a distance equal to its diameter D: 
the mean speed of the large scale roll being l.Q cms~^. This is also in agreement with 
Lagrangian measurements made using density matched tracers: the velocity spectrum was 
found to follow the expected Lagrangian behavior up to a frequency set be the particle's 
size and the flow characteristic velocity scale [sl. These interpretations are consistent with 
qualitative observations using Schlieren visualization, which reveal that temperature plumes 
of size D or larger fully entrain the smart particle, while smaller plumes do so only partially. 

Finally, as shown in figure 3c, the histogram of temperature is non-Gaussian. The non- 
Gaussian character of the PDF is in agreement with Eulerian measurements using fixed 



temperature probes [13|. Other recent experimental studies have also shown that the distri- 
bution of intense temperature plumes in the form of mushrooms has non Gaussian tails ^. 
We note that the asymmetry of the PDF is probably due to imperfect density matching 
between the fluid and the particle, so that it tends to stay slightly longer near the bottom 
plate. 
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C. Heat flux measurement 



With the assumption that the particle may indeed act as a Lagrangian tracer in the 
convective flow, it is natural and interesting to discuss the heat carried by the flow in terms 
of the temperature measured by the particle. To this end, we define an instantaneous 
Lagrangian Nusselt number as: 

Nu\t) = l + -^e'it).v,{t) , (2) 

where 6'{t) = 6{t) — 6 is the particle temperature variation from its time averaged tempera- 
ture 6, and Vz(t) its vertical velocity. Integrating this quantity over all the fluid particles in 
any horizontal plane recovers the traditional Nusselt number. 

The time series of Nu^ is shown in figure 4a. Comparing this time series with the tra- 
jectory in figure 2(c), we see that the heat transfer intermittently spikes to very large values 
often as the particle moves away from the end plates and mixes with fluid at intermediate 
heights in the vessel. The power spectrum in figure 4b shows that the periodicity of the 
large scale convection roll is diminished and strong fluctuations occur over a broader range 
of (slow) scales. The spectrum shows a fast, decrease in the same frequency range that 
the temperature showed /~^. As discussed in the previous section, this scaling is consistent 
with our understanding of velocity spectra (C^ ~ /~^) resulting in Nu^ ~ O^v"^ ~ f~^. 
The PDF of Nu^, figure 4c is skewed toward positive values, as expected because coherent 
convective motion are associated with either hot fluid rising or cold fluid sinking, so that 
in each case 9' ■ > 0. The much less probable events with 9' ■ < correspond, for 
instance, to the rise of the particle when it is colder than its environment. This is unlikely 
but not impossible as the particle may be trapped in the swirls of a turbulent plume (e.g. 
the flow in the "cap" of the mushroom shaped plume is opposite the flow in the "stem"). 
The most probable value of the PDF is zero, indicating the particle spends a large amount 
of time at the mean fluid temperature in the bulk of the cell. The mean value, Nu^, is 
more surprising: at 335.4 it is roughly twice the global, Eulerian, Nusselt value computed 
from power input and velocity differences. Such an increased mean may be due to the fact 
that the particle does not sample the flow uniformly, but may instead get more advection 
from intense thermal plumes. The dramatically non-Gaussian form of the PDF as well as 
scale-by-scale statistics of the local flux requires further study. This analysis will be reported 
elsewhere. 
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IV. CONCLUSIONS 



We report on new instrumentation for making fluid mechanics measurements in tfie ref- 
erence frame of particles passively advected by flow motion. These smart particles show 
promise for a fresh perspective in turbulence research as well as two phase flows. The rapid 
development of MEMS and microfluidic components, which may be incorporated in a smart 
particle, opens a diverse range of potential applications including granular flows and reacting 
flows in chemical and biological systems. 

We demonstrate the capabilities of smart particles in the simple case of Lagrangian 
temperature measurements in thermal convection. We report several original findings. First, 
we observe /"^ power law scaling behavior for the temperature spectrum, which suggests 
the temperature acts as a passive scalar advected by fully developed, Kolmogorov-like, 
turbulence. This result supports the hypothesis that outside the boundary layers, where 



buoyancy sets th^ 
bulk. (e.g. Q 



Q, [l5 



_iuid into motion, the temperature is passively mixed in the turbulent 
161]). Strong temporal intermittency in Lagrangian heat transport is 
revealed in the extremely heavy-tailed probability distribution functions, while the fact 
that the mean Lagrangian heat transport is much higher than the Eulerian average suggests 
strong spatial inhomogeneity as well. Further and deeper investigations of these features 
will be presented soon elsewhere. 
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FIG. 1: (a) Block diagram of measurement system, (b) Photo of the smart particle before closing 
the capsule. 
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FIG. 2: (a) Cartoon of the convection flow. The black disc represents the smart particle. (b,c) 
Trajectory of the particle with temperature (b) and heat flux (c) encoded in the color. 
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FIG. 3: Temperature measurements: (a) time series for the trajectory in figure 2(a). The mean 
temperature is 27.887 degC and the standard deviation is 0.451 degC; (b) corresponding power 
spectrum; (c) Temperature probabihty density function (circles), compared to a Gaussian (sohd 
hne) . 
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FIG. 4: Heat flux measurements: (a) time series of Nu^[t) for the trajectory in figure 2(a) 
corresponding power spectrum; (c) Probability density function. 



14 



